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We propose a fault tolerant loading scheme to produce an array of fermions in an optical lattice 
of the high fidelity required for applications in quantum information processing and the modelling 
of strongly correlated systems. A cold reservoir of Fermions plays a dual role as a source of atoms to 
be loaded into the lattice via a Raman process and as a heat bath for sympathetic cooling of lattice 
atoms. Atoms are initially transferred into an excited motional state in each lattice site, and then 
decay to the motional ground state, creating particle-hole pairs in the reservoir. Atoms transferred 
into the ground motional level are no longer coupled back to the reservoir, and doubly occupied sites 
in the motional ground state are prevented by Pauli blocking. This scheme has strong conceptual 
connections with optical pumping, and can be extended to load high-fidelity patterns of atoms. 

PACS numbers: 03.67.Lx, 42.50.-p, 03.75.Ss 



I. INTRODUCTION 

High-precision control of cold atoms in optical lattices 
has found many potential applications in recent years, 
especially in the implementation of quantum informa- 
tion processing and the modelling of strongly correlated 
condensed matter systems pj. These applications have 
been fuelled by experimental techniques which enable en- 
gineering of lattice models with sensitive control over lat- 
tice parameters 0- 0- QJ • independent control for differ- 
ent internal spin states and control of interactions 
between atoms via Feshbach resonances 0, Q • 

For high precision applications, initial state prepara- 
tion will play a key role in addition to such control of 
Hamiltonian parameters Quantum computing appli- 
cations generally require an initial register with exactly 
one atom per lattice site , and observation of interest- 
ing effects in strongly correlated systems often requires 
the initial spatial patterns of atoms or states with pre- 
cisely chosen filling factors [lOj. 

The first step in preparation of such states is often 
adiabatically increasing the lattice potential, making use 
of repulsive onsite interactions for bosons [TI| or Pauli 
blocking for fermions 01 to load essentially one atom on 
each lattice site. However, experimental imperfections 
will generally lead to non-negligible errors in the result- 
ing states. This can be improved upon by coherently 
filteringa state with a filling factor initially greater than 
one [8lll3l|. or potentially by schemes involving individual 
addressing and precise measurement of the occupation 
in individual lattice sites [fH Il5j . Whilst these meth- 
ods can, in principle, produce high fidelity initial states, 
each of them relies either on the perfect experimental 
implementation of a single-shot coherent process or on 
perfect measurements to avoid defects in the final state. 
In this article we propose a fault-tolerant loading scheme 
in which the state being prepared always improves in 
time. The key idea is the addition of a dissipative element 
to the loading process, in contrast to previous schemes, 



which rely on coherent transfer or perfect measurements. 
As we will see below, this dissipative element plays a sim- 
ilar role in our scheme to that of spontaneous emissions 
in optical pumping. 

Motivated by advances in experiments with cold 
fermions [H El Q 13 IH EJ El El , our scheme is 
designed to produce a regular patterned array of fermions 
in an optical lattice. Fermions have a natural advan- 
tage in initialising atomic qubit registers because Pauli- 
blocking prevents doubly-occupied sites, and most of 
the techniques illustrated using bosons in quantum com- 
puting proposals apply equally to fermions. Fermionic 
species are also of special interest in the simulation of 
condensed matter systems |25j . 
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FIG. 1: Laser-assisted dissipative loading of fermions in an op- 
tical lattice: (a) Atoms are coupled from an external reservoir 
(in internal state \b)) into an excited motional state in the lat- 
tice (internal state |o) via a Raman process; (b) These atoms 
are cooled to the ground motional level via collisional interac- 
tion with the reservoir atoms, and doubly-occupied sites in the 
ground level are prevented by Pauli blocking; (c) Remaining 
atoms in the excited motional levels are removed by carefully 
detuning the Raman coupling above the Fermi Energy. 

The setup for our scheme is illustrated in Fig. ^ 
Atoms in an internal state \b) do not couple to the lattice 
lasers, and form a cold Fermi reservoir, which will play 
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the dual role of a source for atoms to transfer into the lat- 
tice, and a bath for cooling lattice atoms. Atoms in the 
reservoir are coupled into an excited motional level in the 
lattice (in internal state |a)) via a coherent laser-induced 
Raman process (Fig. QJi) [2(j. These atoms are then 
cooled sympathetically by the reservoir atoms via colli- 
sional interactions, and will decay to the motional ground 
state together with creation of a particle-hole pair in the 
reservoir (Fig. ^p). This is analogous to the sympa- 
thetic cooling process previously presented for a bosonic 
reservoir in |27| . Double occupancy in the ground state 
is prevented by Pauli blocking (Fig. QJd), and atoms in 
the ground state are not coupled back to the reservoir 
because the Raman process is far off resonance, so the 
occupation of the lowest motional level always increases 
in time. Additional atoms remaining in excited states at 
the end of the process can then be removed by a care- 
ful adiabatic detuning and switching off of the coupling 
lasers (Fig. ^). 

Such dissipative transfer of atoms into a desired dark 
state is strongly reminiscent of optical pumping, in which 
atoms are excited by a laser, and undergo spontaneous 
emissions into a desired state which does not couple to 
the laser field. The net result of this is to transfer entropy 
from the atomic system into the "reservoir" (the vacuum 
modes of the radiation field) in order to produce a single 
pure electronic state from an initial mixed state. Here 
the creation of an excitation in the reservoir replaces the 
spontaneous emission event, placing the atom in a state 
where it is not coupled by the Raman process, and lead- 
ing to the production of our final pure state, namely a 
high fidelity array of one atom in each lattice site (or a 
pattern of occupied and unoccupied sites). 

We note, in addition, that the purely coherent laser- 
assisted loading could be used as a stand-alone technique 
to load the lattice, and could produce high fidelity states 
if used iteratively, together with cooling of the Fermi 
reservoir. Such cooling would fill holes produced in the 
previous loading step, so that Pauli blocking would pre- 
vent a net transfer of atoms from the lattice to the reser- 
voir, thus ensuring that the filling factor in the lattice is 
improved in each step. 

The detailed analysis of this dissipative loading process 
is divided into two parts. Coherent laser-assisted load- 
ing of atoms into the entire lattice in a single addressed 
motional band is discussed in section |HJ and the dissi- 
pative transfer of atoms to the lowest motional band is 
analysed in section ITTT1 The combination of these two el- 
ements into the overall scheme is then detailed in section 

El 



producing the optical lattice. They form a Fermi gas con- 
taining N atoms with a density «3D = N/V in a volume 

V, with Fermi energy (h = 1) ep = (67t 2 ?13d) 2 ^ 3 /2m, 
where m is the mass of the atoms. The internal state 
| b) is coupled to a different internal state |a), which is 
trapped by a deep three dimensional optical lattice po- 
tential V a {x) via a Raman transition. 

Our goal is to couple atoms b into the lattice and to 
achieve an average occupation of fermions close to one in 
all lattice sites in one chosen motional band, without cou- 
pling to other motional levels. This should be achieved 
in a time, where no atoms are allowed to tunnel between 
different lattice sites and no loss of atoms occurs (due to, 
e.g., spontaneous emission events leading to additional 
internal states etc.). 

A. The Model 

The total Hamiltonian of this system is given by 

H = H a +H b + H RC , (1) 

where the Hamiltonians for the atoms a in the optical 
lattice and for the atoms b forming the reservoir are 

H a = J d'Wl(x) + K(x)) 4(x), (2) 

and 

H b = J d 3 ^(x) &(x), (3) 

respectively, in which the anticommuting field operators 
^|(x) create a fermion in the internal state i £ {a, b} at 
the position x. 




FIG. 2: Reservoir atoms with energy e are resonantly coupled 
to the first excited Bloch band of the lattice via a Raman laser 
with two photon Rabi frequency SI. The resonant energy e 
is experimentally tunable via the Raman detuning, and the 
energy separation of the two Bloch bands is denoted by u>. 



II. LASER- ASSISTED LOADING 

We begin by studying the coupling of the atoms form- 
ing the reservoir into the optical lattice via a Raman pro- 
cess, as shown in Fig.^ The atoms in the Reservoir are in 
an internal state \b), which does not couple to the lasers 



The two internal states are coupled via a Raman pro- 
cess described by the Hamiltonian 



H RC = I d 



2 (4 t ( X )^a(x)+h.C. 

+A^( x )^ a ( x ) 



(4) 
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with the Raman detuning A and the effective (two pho- 
ton) Rabi frequency f2, where we have assumed running 
waves with the same wave vectors for two lasers produc- 
ing the Raman coupling. 

We expand the field operators for the free fermions in 
the reservoir as plane waves and the field operators for 
the lattice atoms in terms of Wannier functions, 



Mx) = J=]Te ikx 6 k , 
v v k 

V^a(x) = 2J W n (x - X Q )a Q , n , 



(5) 



where, creates a reservoir atom with momentum k, 
aj, n is the creation operator for an atom in lattice site 
a and motional state with n = (n x ,n y ,n z ) in the deep 
three dimensional optical lattice, for which u; n (x — x a ), 
denotes the corresponding Wannier function. 
Inserting into Eqs. ©-(0J we obtain 



H, 



ek6 k & k 



k,a,n 

where the single particle energy of a reservoir atom with 
momentum k is e k = |k| 2 /2m and the energy of a lattice 
atom in the motional state n is given by 

w„ = J d 3 a;w n (x) {-—^ + V a (x)J %(x). (7) 

As we are dealing with very deep optical lattices, tunnel- 
ing between different lattices sites is strongly suppressed 
and has thus been neglected. The Raman coupling pa- 
rameter i?k.n can be written as 



1 



d 3 o 



-ikx 



(8) 



For our deep optical lattices without tunneling between 
different sites, the periodic lattice is equivalent to an 
array of independent microtraps, where each individual 
trap is well approximated by a harmonic oscillator. The 
Wannier functions w n (x. — x Q ) can then be approximated 
by harmonic oscillator eigenfunctions </> n (x — x Q ) of the 
n-th oscillator level in lattice site a. This approxima- 
tion allows us to calculate the coupling parameters i? k , n 
from the Fermi reservoir to the optical lattice explicitly. 
For an isotropic three dimensional lattice (where the fre- 
quency uj of each oscillator is given by uj w uj n — ui n ~i 
and n = (n x + n y + n z )) the couplings to the lowest and 
first three (degenerate) excited motional states are given 
by 



Rk.O = —f= v 



3 / 4 (Rn*\ 1/2 p" k2 «o/2 



where ao = yjl/rmv denotes the size of the harmonic 
oscillator ground state, and the index x, y, z labels the 
coupling to the three degenerate states of the first excited 
oscillator level. 

The characteristics of the coherent loading procedure 
strongly depend on the interplay between the (experi- 
mentally adjustable) parameters: the detuning A and 
two photon Rabi frequency Q of the lasers producing the 
Raman coupling, the Fermi energy ep and the separation 
uj of the oscillator levels. The Raman detuning can be 
adjusted to address different states in the Fermi sea and 
different motional states in the lattice. In the following 
we write A = — 5w/2 + e to indicate the resonant coupling 
of reservoir atoms with energy e to the n — 1 motional 
states of each lattice site. We note that it is straight- 
forward to address other motional states in the lattice 
(e.g. to directly load the lowest level) by adjusting the 
detuning A. However, as we will later use the transfer 
of the atoms from the reservoir to the lattice as a first 
step of an indirect loading of the lowest motional states 
as described in the introduction, we choose the transfer 
to the first excited motional state here. To be able to 
selectively fill the first excited oscillator levels, the condi- 
tions (cf — e) <§C uj and e <C u> (and consequently tp "C oj) 
have to be fulfilled, in order to avoid unwanted coupling 
to higher excited and to the lowest motional state, re- 
spectively. 



B. The Fast and Slow Loading Regimes 

The physics of the loading process allows us to iden- 
tify two different loading limits: (1) the "fast loading 
regime" , where 



and (2) the "slow regime" , where 



(10) 



(11) 



Below we will see that our goal to selectively fill a certain 
motional state without coupling to other states can only 
be achieved in the slow loading regime, but to obtain 
more insight into the physics of the loading dynamics it 
is instructive to discuss both regimes. 

In the fast loading regime, where the Rabi frequency £1 
is the largest frequency scale in the system, the loading 
is performed in a very short time T ~ n/tt <C clq/vf, 
with vf = \/2eir/M the Fermi velocity, where atoms in 
the Fermi reservoir do not move significantly during the 
loading on a lengthscale given by the size ao of the har- 
monic oscillator ground state. The Wannier modes in the 
lattice then couple to localized reservoir fermions at each 
site, and thus the dynamics for different sites decouple. 
Given there is at least one fermion in the reservoir per 
size ao of the ground state in each lattice site during the 
loading, i.e., given the density of the reservoir atoms 



(9) 



"3D > 1/ajj, 



(12) 
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each n = 1 motional state in each lattice site can be 
filled with at least one atom from the reservoir by ap- 
plying a 7r-pulse Q.T = n. For an optical lattice with 
ujjl-K ~ 50kHz the required densities of the Fermi gas 
are > 3 x 10 15 cm~ 3 for 40 K and for deeper lattices 
the required densities are even higher. The condition l|12(l 
for the density of the reservoir can be expressed in terms 
of energies as 



(13) 



This inequality violates the condition tj? < w, which is 
necessary to be able to selectively address individual mo- 
tional states. Consequently, unwanted population will be 
transferred to additional motional states in this loading 
limit, which would have to be carefully removed after the 
loading process. 

In the slow loading regime, where condition is 
fulfilled, the atoms in the Fermi reservoir are no longer 
frozen during the loading process, but are allowed to 
move with respect to the lattice during the loading. This 
is now performed in a time T 3> X/2vp, where A/2 is the 
lattice spacing. Consequently, the density condition itT^) 
can be relaxed to 



"-3D 



> 1, 



(14) 



i.e., we only need one atom in the reservoir per lattice 
site to be able to efficiently fill the lattice. For typical 
experimental parameters A ~ 800nm for 40 K this results 
in the condition n^Y) ^ 10 13 cm~ 3 , which has already been 
achieved in current experiments (e.g. |2S[V The density 
condition (|14|) expressed in terms of energies now reads 



( F 



> 



2/3 



(15) 



with ujr — 2ir 2 jm\ 2 the recoil frequency. As u>r <C oj 
for a deep optical lattice, the condition tp <C ui can be 
fulfilled in this loading limit, and as fl <C ui, individual 
motional states in each site can be addressed. In the 
following we will investigate these two extreme limits and 
the intermediate regime in detail. 



C. Analysis of the Loading Regimes 

1. Fast Loading Regime 

In this regime, where the motion of the atoms in the 
reservoir is frozen on the scale during the transfer, the 
physics is essentially an on-site coupling and transfer. We 
thus find it useful to expand the modes in the reservoir 
in terms of localized Wannier functions corresponding to 
the lattice. Such an expansion of the reservoir modes 
arises naturally from the definition of the matrix elements 
-Rk,ne~ lkXa (Eq. @) and allows us to write 



n 



#RC = J^2 ( B <L a a,n + h.C.) 



(16) 



where B a ll = ^ k i?k. n &ke~ Q is the mode correspond- 
ing to the Wannier function w n (x — x a 
collective modes fulfill 



Note that these 



{B a ,n, Bp m } = 5 at p5 n 



(17) 



(where 8 denotes the Kronecker Delta), i.e., modes corre- 
sponding to different lattice sites or to different motional 
states are orthogonal. Furthermore, in the fast regime we 
can neglect the first two terms H a and Hi, in the Hamilto- 
nian due to the condition Eq. I|l()(l during the loading 
time T ~ 7r/J7 and the total Hamiltonian can be approx- 
imated by H sa Hrc- The sites thus decouple, and the 
loading process at each site proceeds independently, but 
with the same Rabi frequency f2 for the coupling. 

We are interested in the time evolution of the ma- 
trix elements of the single particle density matrix, i.e., 
( a a,n a /3,m), (4,n^,m) and (B^B >m ) . In the fast 
loading regime, where H « -ffrtc, the respective ma- 
trix elements can be calculated analytically from the 
Schrodinger equation with the Hamiltonian Eq. 116|l . 
and we find for states, where f2/2 ^> |(n — l)w|, i.e., 
t < T <C a /v F 



( a Q,n a /3,m(£)) = <5a,/3<5n,: 



and 



(Bl n B , m (t)) 



CI 



Sa.^S^m COS — i, 



(18) 



(19) 



for the time evolution of the occupation of the modes 
in the lattice and in the Fermi sea, respectively. These 
expressions assume that the lattice modes are initially 
empty and the corresponding modes in the Fermi sea 
are initially filled. If the Fermi sea is initially filled up 
to €p, then this assumption is fulfilled for any a and n 
for which that each mode B a ll contains contributions 
only from states with energy below ep . Thus, in the 
fast loading regime the occupation in the lowest and first 
excited motional state undergoes Rabi-oscillations at a 
Rabi frequency Cl, and provided the density is sufficiently 
high, the lattice can be efficiently filled by applying a ir- 
pulse, 



(20) 



with loading time T. Atoms will also be coupled to other 
motional states in the lattice with the resulting filling 
factors depending on the density of the reservoir gas and 
the actual value of the Rabi frequency Cl. 

To model the full loading dynamics we use numeri- 
cal simulations of the dynamics generated by the Hamil- 
tonian ijTjl. In these simulations we only consider the 
lowest two motional states for simplicity, but all results 
are easily extended to more motional states. Also, the 
simulations are one dimensional, which means that the 
excited oscillator state with n = 1 is no longer degener- 
ate. Because couplings to motional excitations in differ- 
ent spatial directions are independent, such simulations 
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FIG. 3: Numerical results in the fast loading limit, (a) The 
time evolution of the occupation of the lowest (n — 0) and 
first excited motional (n = 1) levels against time in dimen- 
sionless units. Note that these lines are indistinguishable, 
(b) 1 — J-„(t — tt/Q) after applying a 7r-pulse, for the low- 
est (n = 0, dotted line) and the first excited (n = 1, solid 
line) motional state, as a function of the dimensionless den- 
sity niDao. Parameters used: N=201 particles in the Fermi 
sea, M=5 lattice sites, Q = 17.8ef, A = — 3cj/2 and in (a) 
uj — 0.1 €f, whereas in (b) u> is varied. 



are representative for loading into each of the three 3D 
modes. 

In Fig. we show the results of our numerical 
simulations of the complete system described by the 
Hamiltonian Eq. (JJJ in the fast loading limit. In the 
upper and lower part we plot the fidelity J- m {t) = 
12 a ( a a m a a,m(t))/M, with M the number of lattice sites, 
of the lowest and first excited Bloch band as a function 
of time in dimensionless units ttF- The numerical results 
are in excellent agreement with the analytical calcula- 
tions (Eq. (|18fl ). as we find oscillations of the fidelity in 
both Bloch bands between zero and J- m {t) > 1 — 1CP 4 
occur with a Rabi frequency tt. In Fig. [2b we analyze the 
scaling of the fidelity in the two bands with the dimen- 
sionless density n^aQ (i.e., with niDdo = \f2ep/oj/'K 
in our one dimensional simulations, with n\p> the one 
dimensional density of the reservoir gas). As expected, 
the fidelity after a ir pulse, i.e., J- m (t = 7r/f2) increases 
with the density, and high fidelity states can be achieved 
for large densities nn>ao > 1. In this and all numerical 
simulations below we have checked that the results are 
independent of the quantization volume, which is much 
smaller than in a real experiment, due to the comparably 
small number of particles in the simulations. 

In Fig. 01 we show how the loading dynamics change 
when approaching the intermediate regime from the fast 
limit, i.e. the scaling of the fidelity with the Rabi fre- 
quency O. In Fig. 0t we show the qualitative behaviour 
of the loading dynamics for typical parameters, in Fig. 0b 
the scaling of the fidelity T m (t — tt/O), m = 0, 1 is shown 
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FIG. 4: Numerical results for fl approaching an intermediate 
regime from the fast loading limit, (a) The time evolution 
of the occupation of the lowest (n = 0, dotted line) and first 
excited (n = 1, solid line) motional state against time in di- 
mensionless units, again for N = 101, M — 5 and for a typical 
set of parameters £1 — 0.72eF, A = —3(^/2 and niDao ~ 2. 
(b) The occupation number for the lowest (n = 0, dotted line) 
and first excited (n — 1, solid line) after applying a 7r-pulse 
versus the dimensionless Rabi frequency Q. 

as a function of the Rabi frequency. These numerical 
simulations show that the Bloch bands still cannot be 
individually addressed, and the fidelity becomes worse if 
the Rabi frequency is decreased. 

Thus, our chosen motional state can, in principle, 
be efficiently filled in this regime on sufficiently fast 
timescales. However, the requirements on the density 
are difficult to achieve experimentally, and occupation in 
other motional states cannot be avoided. As a result in 
this regime we obtain no significant advantage over tradi- 
tional loading mechanisms such as adiabatically turning 
on the lattice. In the next section we will investigate the 
slow loading regime. In this limit these problems do not 
exist and we are able to selectively load a single energy 
level efficiently 



2. Slow Loading Regime 

In this regime, transport is significant during the load- 
ing, and the system dynamics are described by the com- 
plete Hamiltonian (JTJ. As the reservoir atoms move be- 
tween Wannier modes during the loading process, it is 
now more convenient to directly use the momentum rep- 
resentation Eq. © to express the coupling Hamiltonian. 

From Eq. © one can see that each lattice site a and 
each motional state n is coupled to many momentum 
modes 6k in the reservoir. However, as £1 <g. ep, effec- 
tively only a subset of momentum modes with energies 
centered around the resonant frequency e = A + 3w/2 is 
coupled to the lattice, whereas the remaining states are 
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far detuned and the transfer is suppressed. The width of 
this effective coupling range depends on both the Rabi 
frequency SI and the matrix elements i?k, n , and an upper 
bound for the width of this range is given by the Rabi 
frequency SI. 

It is convenient to rewrite the coupling Hamiltonian of 
Eq. © as 



RC 



E 



-ikx Q 



h.c. 



(21) 



from which we can see that each momentum mode in 
the reservoir couples to a collective mode J2 a elVk c " a a.n 
in the lattice. To fill the lattice it is necessary that the 
range of states in the reservoir couples to at least M 
orthogonal collective modes in the lattice. Writing the 
phase as 



k x a 
kp~~\/2 



(22) 



we see that it is necessary to couple a range of states with 
width of at least kFyUii/tF in momentum space to the 
lattice to fill M lattice sites. In the slow regime, where 
SI <C ep and furthermore lor < cf (from the density 
condition l|15|l. the recoil frequency will typically exceed 
the Rabi frequency, i.e., SI < lur. As only states within a 
range e ± SI are coupled to the lattice, the lattice cannot 
be filled efficiently for a constant e. 

Thus to achieve a high population in the desired mo- 
tional state of each lattice site we must sweep the reso- 
nant frequency e through a range of at least ujr, scanning 
through many modes. Such a procedure also has the ad- 
vantage that as we only couple to a narrow range in the 
Fermi sea at any one time, the reverse process of trans- 
ferring particles from the lattice to the Fermi sea will be 
suppressed by Pauli blocking. In our numerical simula- 
tions we linearly sweep the detuning from e = to e = ep 
in a loading time T. 

We are interested in the time evolution of the matrix 
elements (a f a n a/3. m ), (a^ n 6 k ), and (fc^iv) of the single 
particle density matrix. For a system described by a 
quadratic Hamiltonian the equations of motion for the 
second order correlation functions can be obtained from 
the (linear) Heisenberg equations (see Appendix 0. As 
the system is described by the quadratic Hamiltonian 
and © , the linear Heisenberg equations for the operators 
a Qi „ and &k have the simple form (again only consider- 
ing the lowest two motional states in a one dimensional 
system) 

. SI ^ # :i_ x , 

a>a,i = -i-^- > y ^ q ,ie Q o q -iea Qj i, 



iek&k, 



.SI 



(23) 



which can be used to efficiently calculate the time evo- 
lution of the desired functions numerically. Note that 
in an isotropic three dimensional lattice again all three 
degenerate n — 1 states will be loaded by sweeping the 
resonant frequency e through the Fermi sea. In practice 
it is also possible to selectively load only a single atom in 
each lattice site by shifting two excited motional states 
out of resonance, choosing an anisotropic lattice with sig- 
nificantly higher oscillator frequencies in two dimensions. 
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FIG. 5: Numerical simulation of the loading dynamics in the 
"slow loading" limit, (a) Occupation of the lowest and first 
excited motional bands as a function of time, showing attain- 
ment of a high fidelity in the excited band. The resonant 
frequency e is swept from ef — > 0, and Q is ramped from 
to 0.45eF, reaching that value at te F — 500. Parameters 
used: N = 81 particles in the Fermi sea, M — 5 lattice sites, 
a) — 5e F and tiidA/2 = 3.4. (b) The final occupation number 
after a loading sweep with constant SI, and e ramped from 
€f — » 0, versus the dimensionless sweep time tfT. Param- 
eters used: N = 81, M = 5, S} — 0.9ef and ui — 10ef and 
nioA/2 = 1.7. 

In Fig. we show numerical results for the time evo- 
lution of the occupation number in the first (upper plot) 
and in the lowest (lower plot) Bloch band as a function 
of time in dimensionless units. Here, SI is slowly switched 
on to reduce the additional holes introduced in the Fermi 
sea by coupling atoms into states above ep. This is an 
example of many possible optimisations to produce high 
filling, and we find the final T\ > 0.99, in a time of the 
order of 10 milliseconds (with u> ~ 2-7T x 100 kHz). In 
Fig. [SJd the occupation of the two motional levels after a 
loading sweep is plotted as a function of the sweep time 
T. These results are not optimised (SI is held constant, 
and we sweep e from tp — > 0), but still produce fideli- 
ties T\ > 0.95 on a timescale of a few milliseconds, and 
we see that the average filling factor increases with the 
loading time. 

It is important to note that whilst high fidelities can be 
obtained by optimising the parameters of the sweep, it 
is not necessary to achieve high filling during this sweep 
in order to produce high fidelities for the overall load- 
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ing scheme. In the full scheme with decay of atoms to 
the ground motional state included, the upper band need 
never be completely filled at any one time, and removal 
of atoms via the decay process will lead to further atoms 
being coupled into the lattice in the upper motional band. 

Due to the condition f2 <C u), unwanted coupling to 
other Bloch bands can be avoided in this regime, by 
choosing e <§; oj (c.f. Fig. [TJ, as the coupling is then 
sufficiently far detuned as demonstrated in the lower two 
plots of Fig. |3J The scaling of the unwanted coupling to 
the lower band is shown in Fig. [SJi, where we plot the 
occupation of the two Bloch bands after a linear sweep 
with tpT = 300 against the ratio ui/tp. 

In Fig. 03 we show the numerical results when ap- 
proaching the intermediate regime, i.e, the scaling of the 
occupation of the two bands after the linear sweep with 
the Rabi frequency. We find that also here high occu- 
pation of the first Bloch band can be achieved, but by 
increasing the Rabi frequency the unwanted coupling to 
the lower band also increases, as can be seen in the lower 
plot of the figure. 
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FIG. 6: In (a) we show the occupation of the lowest and 
first excited Bloch band after a linear loading sweep from 
e — £f —> with epT = 300 versus the band separation ui 
in dimensionless units. Part (b) shows the loading dynamics 
approaching an intermediate regime from the slow regime: 
We plot the occupation of the lowest and first excited Bloch 
band after a linear loading sweep, with efT = 300 against 
Q. Parameters used: N=81 particles in the Fermi Sea, M=5 
lattice sites and nioA/2 = 1.7, in (a) £1 = 0.9eF and in (b) 

In summary, in the "slow loading" regime high fidelity 
loading of the n = 1 motional level can be achieved on 
timescales much shorter than those on which atoms are 
lost from the lattice by sweeping the resonant coupling 
frequency e through the Fermi sea. This loading mecha- 
nism gives us the significant advantage over simple load- 
ing techniques such as adiabatically increasing the lattice 
depth that we can address a particular energy level in the 
lattice, whilst not coupling to levels that are sufficiently 
far detuned. This property can also be used to load pat- 



terns of atoms, because if a superlattice is applied, then 
the energy of certain lattice sites can be shifted out of 
resonance with the Raman process, so that no atoms are 
coupled into these sites. 

In the next section we will discuss the cooling of atoms 
in higher motional levels to the ground state, which re- 
moves the atoms from the motional state being cou- 
pled from the reservoir. Together with Pauli blocking 
of modes in the lattice, this allows us to make the overall 
loading process fault-tolerant. As an additional remark, 
though, we note that this laser-assisted loading of a se- 
lected energy level in the lattice could be used as a stand- 
alone technique to load the lattice, e.g., coupling atoms 
directly into the ground motional state. (In order to load 
an excited motional state in this manner, interaction of 
atoms in the lattice and atoms in the reservoir must be 
made very small on the timescale of the loading process, 
e.g., by using a Feshbach resonance, in order to avoid 
decay of the atoms into the ground state). This process 
on its own is not as robust as the procedure we obtain 
by including a dissipative element in the loading scheme, 
which will be discussed in the next section. However, 
reasonably high fidelities could still be obtained with this 
method alone, especially if the method was applied itera- 
tively, cooling the Fermi reservoir between each two steps. 
Net transfer of atoms already in the lattice back to the 
reservoir would be prevented in each step by Pauli block- 
ing in the filled Fermi sea. Note again that as with the full 
dissipative loading scheme, a single sweep would also not 
need to completely fill the upper band. The dissipative 
element discussed in the next section allows for the pro- 
duction of an arbitrarily high-fidelity state without the 
requirement of iteratively cooling the Fermi reservoir. 



III. DISSIPATIVE TRANSFER: COOLING 
ATOMS TO THE LOWEST BAND 

The second stage of the loading process is cooling 
atoms in an excited motional state to the ground state 
via interaction with the reservoir gas. This is closely re- 
lated to the cooling process with a bosonic reservoir in 
|27j . The external gas here plays the role of an effec- 
tive T — heat bath for the lattice atoms, and ground 
state cooling is achieved on timescales much shorter than 
atoms are lost from the lattice. 

We consider the coupling of lattice atoms a via a colli- 
sional interaction to the atoms b in the reservoir so that 
the system is described by the Hamiltonian 



H — H„ + Hi, + H\, 



(24) 



where the collisional interaction, H{ nt , between two 
fermions is the usual density-density interaction 



with 



ffmt =9 I d^^(x)Va(x)^ T (x)^(x), (25) 
: 4:Tra s /m and a s the s-wave scattering length. 
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Expanding the field operators as described in the previ- 
ous section we obtain 



with 



-Hmt = 



E 5a;n,n' 6 k 5 k'4, n aa,n' 



(26) 



k,k' 



which is local in each lattice site because of the small 
overlap between Wannier functions for neighboring sites 
in a deep lattice, with 



k,k' 
?a,n,n' 



v e 



,(k'-k) 



d 3 xe ix ( k'-k; 



Wn(x)Wn' (x). 



(27) 



k,k' 



Each g a ' a n , describes a scattering process in which a 
particle-hole pair is created in the reservoir by scattering 
an atom from momentum state k — > k', combined with 
the transition of an atom at site a from motional state 
n — * n'. 

If the transition in the lattice is from a higher en- 
ergy mode to a lower energy mode, this corresponds to 
a cooling transition, whereas the reverse process consti- 
tutes heating. As the initial temperature of the reservoir 
ksT C £f < w, the heating processes will be, at least 
initially, insignificant, as few reservoir atoms will exist 
with sufficient energy to excite an atom in the lattice. If 
the number of atoms in the reservoir is large compared 
to the number of sites in the lattice (N 3> M), then the 
rate of heating processes due to interaction with previ- 
ously excited atoms will be small compared to cooling 
processes due to interaction with atoms remaining be- 
low the Fermi energy ep. Because the cooling processes 
in different lattice sites couple to different modes, and 
therefore are incoherent, the reservoir can then be treated 
throughout the process approximately as a T = bath. 

This can be further enhanced in two ways. Firstly, 
in an experiment in which the reservoir gas is confined 
in a weak harmonic trap, particles with sufficiently large 
energies can be allowed to escape from the trap. The 
large separation of the Bloch band w, and correspond- 
ing excitation energy will then cause many excited reser- 
voir atoms to leave the trap, providing effective evapo- 
rative cooling during the process. Secondly, the lattice 
depth could be modulated during the experiment, so that 
the excitation energy changes, decreasing the probabil- 
ity that atoms are heated by previously excited reservoir 
atoms. 

The cooling dynamics are then described in the Born- 
Markov approximation by a Master equation for the re- 
duced density operator p for the atoms in the lattice. If 
we consider coupling of atoms from the first excited mo- 
tional levels n e {(1, 0, 0), (0, 1, 0), (0, 0, 1)} to the ground 
state, the resulting master equation (derived in Appendix 
□3 is 

9= E ^( 2 ^,nMkn-4,n^,nP-p< n ^,n 



a,/3.n 



k,k' 
k' >k 



g 2 n 3 ym 2 
y/2 3e 



7ra 



(29) 



Here, the jump operator A a jl — u a u a , n 
cooling of a lattice atom in site a from the first excited 
motional level n to the ground state. These results are 
obtained by calculating the integral over momenta in the 
Fermi sea to lowest order in cf/uj. 

The approximation in the second line of Eq. i|29|) . in 
which neglect off diagonal terms a / /3 amounts to the 
approximation that the coherence length of the Fermi 
reservoir is much shorter than the lattice spacing. This is 
true provided that the wavelength of the emitted particle 
excitation, ^2tt 2 / (moj), is much shorter than the lattice 
spacing, i.e., ojr/oj <C 1. This is consistent with the 
previous approximation that the lattice is so deep that we 
can neglect tunnelling between neighbouring sites. This 
can be seen directly when these off-diagonal terms are 
calculated, as for large lor/lo they decay (to lowest order 
in ep -C u>) as 



sm(iTy/Lj/u R \a - fi\) 
ir^u)/u R \a - f3\ 



(30) 



(28) 



This effect is analogous to the spontaneous emission of 
two excited atoms which are separated spatially by more 
than one wavelength of the photons they emit. In this 
case, the atoms can be treated as coupling to two inde- 
pendent reservoirs, and effects of super- and sub-radiance 
do not play a role. 

For typical experimental values n^o ~ 10 14 cm~ 3 and 
a s = HAaB, for 40 K as given in [2^, with the Bohr radius 
as and a deep optical lattice with lu/2tt ~ 100kHz, we 
find a decay rate r/27r ~ 3.6kHz. Thus, cooling can again 
be achieved fast enough, as this rate is much faster than 
typical loss rates of the lattice atoms. Note, that this 
value of the decay rate can be made even larger e.g. by 
tuning the scattering length a s via a Feshbach resonance, 
as r cx a s , by increasing the density of the external gas 
or by increasing the lattice depth. 

In summary we have shown that for a cold reservoir 
gas with sufficiently many atoms fast ground state cool- 
ing of lattice atoms can be achieved with the dissipative 
coupling of the lattice to the reservoir. The necessary 
experimental parameters have already been achieved in 
real experiments, and the cooling rates are tunable via 
the scattering length and the density of the reservoir gas. 



IV. COMBINED PROCESS 

The combination of the cooling process with laser- 
assisted loading in the limit Q <C to,ep will give a final 
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high-fidelity state in the lowest motional level. The pri- 
mary role of the dissipative element is to transfer atoms 
into a state in which they are not coupled back to the 
Fermi reservoir, which is made possible because of the 
selective addressing of the n = 1 motional levels in this 
regime. Multiple occupation of a single site in the lowest 
motional state is forbidden due to Pauli-blocking, and 
thus the lowest motional state is monotonically filled; 
the filling factor and hence the fidelity of the state be- 
ing prepared always improving in time. Again, patterns 
of atoms may be loaded in the lowest state by using a 
superlattice to shift the energy of the n = 1 motional 
level out of resonance with the Raman process in partic- 
ular sites, preventing atoms from being coupled from the 
Fermi reservoir into those sites. This energy shift will also 
further suppress tunnelling of atoms from neighbouring 
sites. 

If the laser-assisted loading and the cooling are car- 
ried out separately, each being performed after the other 
in iterative steps, then from the analysis of sections ITT1 
and 11111 we see that an arbitrarily high fidelity final state 
can be obtained. This pulsed scheme gives us an upper 
bound on the timescale for loading a state of given fi- 
delity, which corresponds to the combination of the two 
individual timescales for laser-assisted loading and cool- 
ing. Provided that the number of atoms in the reservoir 
is much larger than the number of lattice sites to be filled 
(N 3> M), and the Markov approximation made in de- 
scribing the cooling dynamics is valid, then there will be 
no adverse effects arising from the loading and cooling 
processes sharing the same reservoir. Thus, we can com- 
bine the two processes into a continuous scheme, which 
in practice will proceed much faster, as the continuous 
evacuation of the excited band due to cooling will also 
speed up the loading process. 

At the end of the loading process we must still ensure 
that the finite occupation of the excited motional levels 
is properly removed. This can be achieved by detuning 
the resonant frequency for the Raman coupling above 
the Fermi energy after the loading sweep, coupling the 
remaining atoms to empty states above the Fermi sea, 
and then switching off the coupling adiabatically. 

The dynamics of the pulsed process are already well 
understood from the analysis of sections |TT1 and ITTT1 To 
illustrate the dynamics of the combined continuous pro- 
cess, we again perform numerical simulations, in which 
we compute the matrix elements of the reduced system 
density operator. The dynamics of the total system in- 
cluding both the laser coupling and the collisional inter- 
action between the optical lattice and the Fermi reservoir 
are described by the full Hamiltonian 

H = H a + Hb + Hue + Hi n t, (31) 

and in the Markov approximation with respect to the 
cooling process, the matrix elements of the system den- 
sity operator can now be calculated from the Master 
equation (|28|) as shown in Appendix IU1 In order to ob- 
tain a closed set of differential equations which can be 



integrated numerically, we use an approximation based 
on Wick's theorem to factorize fourth order correlation 
functions into second order correlation functions (see ap- 
pendix 
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FIG. 7: The occupation of (a) the lowest (n = 0) and (b) first 
excited (n — 1) motional level for the continuous combination 
of coherent loading in the slow regime and dissipative cooling. 
The resonant frequency is swept from e = to e = 4e F and 
the Raman coupling is switched off adiabatically. Parameters 
used: N = 101, M = 5, fi/2 = 0.45e F , n m \/2 = 1.4, u = 
10e F and T = 0.1e F . 

In Fig. we plot the time evolution of the occupation 
of the two motional levels in the continuous regime as a 
function of time. In Fig. [7]} we see that we indeed achieve 
a high occupation of the lowest motional level from the 
combination of laser-assisted coupling to the excited mo- 
tional level in the regime Q <C oj,6f and cooling to the 
ground state. For the typical values given in the figure 
caption, the loading time T for a state with ~ 1 — 10~ 4 
is again on the order of a few milliseconds. This required 
loading time can be further decreased by tuning T via the 
density of the external gas or the strength of the colli- 
sional interaction between atoms in the lattice and atoms 
in the reservoir. 

From Fig. [7jp we see that as we fill the lower motional 
level, the filling in the upper level is depleted, and as we 
continue to tune the lasers so that this level is coupled 
to states in the reservoir above the Fermi energy sp, the 
remaining population in this level is removed. 

As a final remark we note that such a procedure could, 
in principle, also be applied to bosons. However, with- 
out Pauli blocking to prevent double-occupation of the 
ground motional level, we rely on the onsite collisional 
shift U to make the Raman coupling of an atom from 
the reservoir into an excited state off-resonant if an atom 
already exists in the ground motional state. Second or- 
der processes occurring at a rate ~ QT/U can still create 
double occupation, so we require Vl, T <C U, and the ad- 
vantage of true fault tolerance is not present as it is for 
fermions. 
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V. SUMMARY 

In conclusion, we have shown that the combination of 
laser-assisted loading of atoms into an excited motional 
state and the cooling of atoms from this motional state 
to the ground level gives a fault-tolerant loading scheme 
to produce high fidelity registers of fermions in an opti- 
cal lattices with one atom per lattice site. Application of 
a superlattice allows this to be extended to generalised 
patterns of atoms, and all of these processes can be com- 
pleted on timescales much faster than those on which 
atoms can be lost from the lattice. The advantage of this 
scheme is that the dissipative transitions in the lattice, 
similar to optical pumping, gives us a process in which 
the fidelity of the final state (in the lowest motional level) 
improves monotonically in time. 
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APPENDIX A: DERIVATION OF THE 
HEISENBERG EQUATIONS FOR COHERENT 
LOADING 

Consider a system, which is described by a 
Hamiltonian quadratic in a set of operators = 
(Oi, 02, ■ ■ ■ , Od). Then the Heisenberg equations of mo- 
tion can be written as 



0(t) = MO(t), 



(Al) 



with a matrix M, and formal solution 0(t) — UO(0) 
with U — exp(Mt). By choosing the initial conditions 
Qj'(O) = 3j,ot we can construct the full time evolution 
matrix U(t) by solving Eqs. IjAlfl . as 

u iia (t) = UtAWM = °<(*)- ( A2 ) 



The time evolution of the second order correlation func- 
tions is then easily calculated as 



£W)tfA/j(*)0to»(o) 



(A3) 



APPENDIX B: DERIVATION OF THE MASTER 
EQUATION 

In the interaction picture, and after making the Born- 
Markov approximation, the master equation for the re- 
duced density operator p of a system which interacts with 
a heat bath via an interaction Hamiltonian H ln t can be 
written as (see e.g. 30] ) 



Kt) =~ J rfTTr B{ [Hint(t), [H int (t 



-).p(t)<8>0B]]}. 

(Bl) 



Here, ps is the bath density operator, and Tre denotes 
the trace over the bath, which is represented by the cold 
Fermi reservoir in our setup. The interaction between the 
Fermi reservoir and the optical lattice system is given by 
the Hamiltonian (|2(j|) , and in the interaction picture with 
respect to the internal dynamics in the lattice and in the 
Fermi reservoir, 



flint (*) 



/ j ■' a . ii . ii 

k.k' 



(B2) 



As the number of atoms in the reservoir exceeds the 
number of lattice sites, N ^S> M, and as in addition the 
the bath has temperature T ~ 0, the reservoir will ap- 
proximately remain in its ground state, i.e., the filled 
Fermi sea throughout the cooling process, and the bath 
correlation functions are approximately given by 

(bl^K^K) ~ 5 kx,k'/k 2 ,k^ + 5ki,k^ki,k a , (B3) 

where ( • ) = Ti b{ • Pb}- 

For t much larger than the correlation time in the bath 
we can let the upper limit of the integral in Eq. (|Blll go 
to oo, and writing / °° +»>{n-n'))T ^ 5 ^ k _ £k , + 

uj(n — n )), we find 



P = 



V 2 

a,/3,n 



(B4) 



with the jump operator A a n — a a ^ n , 

T a ,p, n =2tt 5a: k ,off k ; k ,o S (uj - e k + e k ,) , (B5) 



k,k ; 

k' >k 



k,k' k,k' 



and where we note that £ k k , So/n.oSaX.o = for 
n 7^ n'. The rate T a ^ tXl rapidly decays with \a — /3|, 
and for each of the three degenerate excited states n € 
{(1, 0, 0), (0, 1, 0), (0, 0, 1)}, the slowest rate of this decay 
is found in the direction of n. In the harmonic oscillator 
approximation we find (for the direction with the slowest 
decay) 



(B6) 
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to first order in e^/w, with the function 

For large £ this result simplifies to the sine function in 
Eq. (|3U|) . For a deep optical lattice where lo ^> ojr, 

F '( 7r \J~~^\ a ~ @\) ~ ^q./3; an d we end up with a standard 
quantum optical master equation (see e.g. |3Qjl. describ- 
ing the decay of an excited lattice atom from each of the 
three degenerate n — 1 states to the n = level at a rate 

r. 



APPENDIX C: EQUATIONS OF MOTION FOR 
COMBINED DYNAMICS 

The time evolution of the expectation value of an ar- 
bitrary system operator O can be calculated from the 

I 



3t(4,i a /3,o) = 

ar(<A> = 

Mv> = 

A closed set of equations can be obtained from Eqs. (|C2|) 
by using Wick's theorem to factorize fourth order corre- 
lation functions into products of second order correlation 
functions according to 

(C1C2C4C4) = (C1C 2 )(C3C 4 ) - (ciC 3 ) (C2C4) + (cic 4 ) (C2C3) , 
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